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Extreme events, a type of collective behavior in complex networked dynamical systems, often can have 
catastrophic consequences. To develop effective strategies to control extreme events is of fundamental 
importance and practical interest. Utilizing transportation dynamics on complex networks as a prototypical 
setting, we find that making the network "mobile" can effectively suppress extreme events. A striking, 
resonance-like phenomenon is uncovered, where an optimal degree of mobility exists for which the 
probability of extreme events is minimized. We derive an analytic theory to understand the mechanism of 
control at a detailed and quantitative level, and validate the theory numerically. Implications of our finding 
to current areas such as cybersecurity are discussed. 

Extreme events occur on a large variety of complex networked systems, examples of which include sudden 
bursts of packet flows in the Internet, jamming in computer or transportation networks, abrupt rise of web 
requests, and power blackouts in electrical power grid, etc. The intrinsic dynamics responsible for extreme 
events are network traffic or flow dynamics'" and the corresponding flux- fluctuation behaviors" Especially, 
abnormally large fluctuations in the flow exceeding the nodal or link capacities can lead to extreme events, which 
either can be triggered by external disturbances or are a type of intrinsically emergent behavior. Extreme events 
can have catastrophic consequences, demanding the articulation and development of effective control strat- 
egies^"''''^''. This problem, despite its uttermost importance to many disciplines, has not been addressed. The 
purposes of this report are to present an efficient and physicaOy implementable method and to derive an analytic 
theory to understand its working. 

This work aims to develop a physically practical control scheme against extreme events emerged from the 
intrinsic dynamic factors by using a simple but extensively representative model that captures the essential 
mechanism of many complex networked systems with information, material, or energy transportation. Within 
this model, as presented in the Methods Section, the random variables essential to the emergence extreme events 
are sums of random variables distributed in a binomial fashion (thin-taUed) which have well-defined first and 
second order moments. Thus, these sums follow the Le'vy-alpha-stable distribution g^ix) with a = 2 ^^'^'', and we 
do not expect the model to characterize all properties observed in reality, such as the "fat-tailed" distributions 
corresponding to a < 2. 

Our principal idea of controlling extreme events is to make the whole network system time-dependent in that 
every node is mobUe and consequently the associated set of links is time-varying. The mobility of the nodes can be 
characterized by a velocity-like parameter. The system can thus be modeled as a mobile network in which agents 
move in random directions but with constant velocity. Between any pair of agents, a link is established to enable 
exchange of packets only when their physical distance is within a pre-defined communication radius. We find 
that, the number of extreme events depends on the velocity and, strikingly, it exhibits a bell-shape, resonant-type 
of functional relation, as exemplified in Fig. 1. That is, there exists an optimal velocity for which the number of 
extreme events can be significantly reduced as compared with that in a static network. This finding suggests a 
practical control scheme: the occurrence of extreme events can be suppressed by making the network mobile with 
proper agent velocity. To understand this striking phenomenon, we develop a detailed analytic theory. To extend 
the control, we find a generalized mobility capable of dramatically lowering the control cost. We also develop a 
low-cost "intermittent strategy" that harnesses extreme events not only in mobile networks, but also in the 
classical ER-random networks^'. Since extreme events encompass catastrophic situations such as large-scale 
breakdown induced by random failures of nodes or attacks, the general principle that dynamic networks can 
be immune to extreme events can have significant applications in the field of security and robustness of complex 
systems. 
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Figure 1 | The number of extreme events versus velocity and mobility. For a network of N = 1200 nodes in a square domain of size L = 10 (arbitrary 
unit), (a), The number of extreme events Hex versus the agent velocity v for W = SON, 80JV, and lOON; (b), fi^x of extreme events versus velocity v when 
nodes execute random walk (R) or deterministic motion (D) within square, circular, and stadium-shaped domains. The areas of the circle and the 
stadium are chosen to be equal to the area of the square LX L. (c). Robustness of control with respect to heterogeneity in the nodal communication range: 
n„ versus vfor /? = — 1, 0, and 1, where /? is a parameter characterizing the distribution of nodal communication range, (d) and (e), robustness of control 
with respect to packet-transportation protocol: n„ versus velocity v when packet generation and annihilation are taken into account for two situations 
where (d) the total number of packets is fixed at Wand (e) the number of newly generated packets at each time step is fixed to be WIN. (f). Generalization 
of mobility, n^^ versus mobility, defined as the probability that a node moves with velocity v, for a number of v values. All simulation results are obtained 
using 100 realizations in T = 1000 time steps. The error bars in from different realizations have relatively small and similar magnitude in (a-e), and it 
monotonously increases with velocity v, as shown in panel (a) [the error bars in (b-e) are omitted for clear visualization]. Here, the definition of an 
extreme event on a node is that its number of packets is at least four standard-deviations above the average. 



We emphasize that our work is focused on controlling extreme 
events caused by the intrinsic randomness of the system, rather than 
external attacks. More sophisticated mechanisms such as failures and 
load redistribution, e.g., as triggered by external attacks, would be an 
interesting topic of study. Instead of analyzing the extreme value 
distributions above the threshold as determined by the extreme value 
theory^*'^'* under the peak-over-threshold setting, our work concen- 
trates on calculating the total number of events that exceed a certain 
threshold for the purpose of control. 

Results 

We consider a square domain of size L X L in which N agents move at 
speed V but in random directions. The communication radius of each 
agent is a<CL, and the total number of packets on the network at any 
given time is W. At time t, the number of packets carried by a node is 
w{t), and each packet is delivered to a randomly selected neighboring 
node at time t + 1. In contrast to previous works on synchronization 
in mobile-agent networks'"' where periodic boundary conditions 
were often used, we use the "hard-waU" boundary condition for the 
consideration that a cross-boundary link connecting two nodes close 
to the same boundary but physically distant from each other is not 
physical in an actual traffic-flow network, which is discussed in a 
later section. We will see that boundary effect plays a crucial role in 
the dynamics of extreme events, a fact that wUl be fuUy exploited in 
our development of the analytic theory. 

The existence of optimal mobility and its robustness. Our exten- 
sive computations indicate that the existence of an optimal moving 
velocity to suppress extreme events is general: it holds under a variety 



of boundary shapes, motion schemes, and even heterogeneously 
distributed communication radii, as shown in Figs. la,b, and c. 
This indicates that the characteristic, non-monotonous relation 
between and velocity v that underlies the control of extreme 
events is in fact robust with respect to variants in the system 
setting such as the geometric shape of the domain, the manner by 
which the individual nodes move, heterogeneity in the nodal 
communication range, and packet transportation protocol. 

As can be seen from Fig. la, different W values do not affect the 
existence of the optimal moving velocity. Figure lb shows, for a 
number of variants in the geometric shape of the domain (square, 
circle, and stadium), the relation between n^-^ and v. We observe that 
the non-monotonous behavior holds for all geometric shapes con- 
sidered. As shown, the particular manner by which nodes in the 
network move in the space has little effect on the non-monotonous 
behavior. In particular, both random walk and deterministic motion 
are considered, where for the latter, random moving directions are 
chosen initially but at time a node moves along the direction of the 
velocity at time t — I. Figure lb indicates that the existence of an 
optimal velocity to minimize the number of extreme events is highly 
robust regardless of the domain geometry, random or deterministic 
movements. 

In realistic systems, mobile agents' ranges are not homogeneous. 
To model heterogeneity in the nodal communication range, we ran- 
domly assign each node a weight jx according to the distribution P(^) 
~ ir^. The area of the nodal communication circle is given by 



A,=NAo- 



EN I: 
1=1 l^l 



where Ao = no' is the average communication 
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Figure 2 | Geometrical setting and the theory result, (a), Illustration of various geometrical quantities used in the analysis: communication radius a, 
distance b from the middle of the square to the target node, distance b„ from the middle of the square to a neighbor of the target node, the width Ab„ of the 
bar region, and the cutoff region of the communication circle, (b). Results of Eq. 2 for db = 0.01 and Ah^ = 0.002 (the main features do not depend 
on the specific values of db and Ah„). 



area over aU nodes, /(,- and jii are the assigned weight values corres- 
ponding to nodes ; and /, and P is an adjustable parameter character- 
izing the heterogeneity of A,. The distribution of A, is thus 

P(Aj) ~ A, For /? > 0, majority of nodes have relatively small 
communication circles except for a few, and the opposite situation 
occurs for P < 0. For /i = 0, the communication circles are all 
identical: namely A,- = Aq and a,- = a for ; e [1, 2, N]. Since 
the degree of a node is proportional to its communication area, it is 
equivalent to a power-law degree distribution of the mobile network 
with the fixed average degree (k) = N/L^Ao- The physical meanings of 
positive and negative values of P are the following. In a realistic 
wireless sensor network, for example, there are usually anchor nodes 
that have much larger communication ranges and play the role of 
central controllers of the whole network. This situation can be char- 
acterized by some positive value of p. There can also be a few nodes 
that malfunction or run low on power supply. These nodes would 
have nearly zero communication ranges while most nodes in the 
network have finite communication ranges. This corresponds to 
the case of a negative P value. Figure Ic shows that the non- 
monotonous behavior in the Wex-v relation holds regardless of 
whether the value of P is positive or negative, indicating the robust- 
ness of mobility-based strategy to harness extreme events in the 
presence of nodal heterogeneity. 

We also consider realistic variants of the packet-transportation 
model by randomly assigning each packet a destination node where 
the packet annihilates upon its arrival at the destination, and new 
packets are generated to balance the annihilation effect. Two packet 
generation schemes are studied: (1) the total amount of packets Wis 
kept constant by generating the annihilated number of packets at 
randomly selected nodes; (2) W/N newly generated packets are ran- 
domly distributed to nodes to keep W statistically constant but with 
fluctuations (the average number of annihilated packets is W/N for 
high velocity). As shown in Figs. Id and e, nearly identical behavior 
in the «ex"i' relation is observed, regardless of the specific packet- 
transportation protocol, providing further support for the robustness 
of our control method. 

While so far we have focused on dynamic networks where agents 
move with certain velocity, the idea of making a network mobile to 
control extreme events can be generalized to different forms of 
"mobility." Here we discuss one variant, in which the mobility is 
defined as the probability that any node moves with velocity v. For 
relatively high values of v, there exists an optimal mobility value at 
which the number of extreme events is minimized, as shown in 
Fig. If Notice that the optimal mobility is near zero, indicating that 
only a small fraction of the nodes actually move and vast majority of 
the nodes remain stationary. The implication is that extreme events 



can be effectively mitigated by moving a few nodes, which is advant- 
ageous from the consideration of control cost. 

Explanation and boundary effect. In the following, we focus on the 
development of analytical understanding of the phenomenon. To 
begin, we note that, since each node receives packets at time t from 
its neighbors that equally distribute all the packets obtained at f — 1 
to all their neighbors, a node's flow w{i) at time t is contributed to by 
its neighbors' flow at time t — 1, denoted by w„(t — 1). In a static 
network, the ensemble averaged flow through a node at last time step 
is proportional to its degree {w{t — 1)) = Wk/{2E). However, in a 
time-varying network, this relation no longer holds. In particular, the 
flow would be independent of the degree under high velocity, leading 
to {w{t — 1)) = W/N. We thus hypothesize the foUowing general 
relation between {w(t — 1)) and node degree for different regimes of 
the velocity v: 

(w(t-l)) = c(v)-fc(0 + <?(v), (1) 

where ^{v) and )}(v) are parameters that depend on the velocity only, 
and the first and second terms correspond to the velocity- dependent 
and velocity-independent flow components, respectively. Generally, 
the quantity c(v) decreases with v, since higher mobility tends to 
yield a weaker flow-degree correlation. However, rjiv) tends to 
increase with v due to the increasing degree of randomness in the 
dynamics. Thus, for v ^ 0, the steady network exhibits a high flow- 
degree correlation and litfle randomness, so we have c(v) = W/{2E) 
and i]{v) = 0. For v ^ =o, the correlation is lost, leading to ^(v) = 0 
and ;j(v) = W/N. Equation (1) thus represents a universal empirical 
law, which we have verified computationally through various models 
(see the Methods section). Conservation of the total number of data 
packets imposes the following constraint on the parameters c(v) and 

niv): W = N^~^^^(w(t-l))P(fc)=>^(v) = [W-N-;7(v)]/(2£). 

Setting l^(v) = 0 (or f;(v) = W/N) yields C(v) = W/{2E) (or ^(v) 
= 0). We see that Eq. (1) contains the results for v = 0 and v ^ °o as 
two limiting cases. 

With the help of Eq. (1), the probabflity for the node at location 
b (illustrated in Fig. 2a) to carry w packet, W{w{t))i„ is obtained. 
This distribution Wi,{w{t)) serves as a centro-symmetric probabil- 
ity field for each dimension of the two-dimensional space. For a 
static node located at bg, the PDF of its flow is H'i,„{w{t)). For any 
moving node (v > 0) initially at bg, the PDF of its flow, ^^{w{t)), 
is contributed by the 4'i,(vf) fields at all the locations it has visited, 
and thus can be expressed as an integral of ^i,{w). In turn, the 
total number of extreme events occurred within T time steps is 
obtained as an integral of 4',3(w(t)) over all possible starting 
location. 
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where (j,3 is the extreme event threshold. Since f;(v) monotonously 
increases with v. Figure 2b shows that Eq. (2) successfully explains 
the numerically observed non-monotonous behavior as exempli- 
fied in Fig. la, providing the theoretical foundation for mobility- 
based control of extreme events. 

Our theory provides a detailed understanding of the effect of velo- 
city on the occurrence of extreme events. To further demonstrate 
this, we exploit the effect of domain boundary on the emergence of 
extreme events by examining the behavior of (w(())^ versus fc, as 
exemplified in Figs. 3al,bl,cl, where a good agreement between 
our theory and simulation for various velocity values can be seen. 

Note that fcj ^ 

radius away from the boundary. For v 
[0, bf\ and decreases with b in the interval 



tf— — —a specifies the position of one communication 

0, {w{t))i, is constant for b e 
II 

bf,- . As V is increased, a 
L 2j 

peak in {w{t))i, arises for b = bp but on the left and right sides of the 
peak, {w{t))y is approximately constant and a decreasing function of 

fc, respectively. For convenience, we call bf. 



the boundary region. 



As a node moves closer to the boundary, more portion of its 
communication circle is cut off, leading to lower values of {k{t)) 
and (w(O)b- For v ^ 0, ^{s{w{t)) is highly localized so that 



*P23(vi'(t)) — >4'(,(iv(t)). For the static network, extreme events 
are approximately equally probable in the domain, except in 
the boundary region, as shown in Fig. 3a2. As velocity is 
increased from zero, 4'fg(w(t)) depends more strongly on 
^fc(M'(0)'s and the peak in (vc(O)fe about bf generates a higher 
value of q<s in the neighborhood region of bp reducing the prob- 
ability of extreme events. In the part of the boundary region 
where low values of (vf(O)b meets with high value of q^, there 
are no extreme events, as shown in Fig. 3b2. For high velocity, 
each node appears in the domain with approximately equal 
probability and ^<s{w(t)) is completely delocalized, so we have 

fi/2 



b-^ h{w[t)) db. In this case, there is little 



difference in the dynamical behavior of nodes. When a node is 
near bp the sudden increase in the flow can cause an extreme 
event, which is the mechanism responsible for high concentra- 
tion of extreme events about b = bp as shown in Fig. 3 c2. 

Figs. 3a3, b3, c3, and a4, b4, c4 show the patterns of the locations of 
the extreme events occurring on nodes for v = 0, 0.5, and 3 within the 
circular and stadium-shaped domains, respectively, where the 
motions of nodes are completely random. For deterministic motions 
of the nodes, the patterns are similar. Regardless of the geometrical 
shape of the domain, a frame of extreme events (FEEs) emerges about 
one communication radius away from the domain boundary for high 
velocity, providing visual support for the basic idea underlying our 
theory: the boundary plays an important role in the emergence of 
extreme events. 
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Figure 3 | Boundary effect, (al), (bl), (cl). Behavior of (w(t))i, versus fefrom direct simulation (blue squares) and predicted by Eq. (9) (solid black lines) 
for V = 0, 0.5 and 3.0, respectively. Note ttiat v = 0.5 is the optimal value of the velocity that minimizes n^^ in Fig. 1 , and v = 3.0 represents a relatively high 
velocity value for the system of size i = 10 and communication radius a= \. For (al), we have )/(v = 0) = 0. Since Eq. (9) indicates a smooth dependence 
of {w(t))y on b and since a detailed interdependence between J7(v) and v is unknown, we find r\(v) = 15 best fits the case in (bl). For (cl), we have 
Jj(v = 3.0) r}(v^ =o) = 50. The red dashed line indicates the position of the peak in the curve of (w(t))6 versus h. (a2), (b2), (c2). Positions of aU extreme 
events occurred in square domain for v = 0, 0.5, and 3.0, respectively. (a3), (b3), (c3). Positions of all extreme events occurred in circular domain 
for V = 0, 0.5, and 3.0, respectively. (a4), (b4), (c4). Positions of aU extreme events occurred in stadium domain for v = 0, 0.5, and 3.0, respectively. AU 
results are obtained from 20 independent statistical realizations. 
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Figure 4 | Effect of long range links, (a), effect of periodic boundary condition (PBC): n^x versus vwhenO% (green circles, no PBC), 30% (aqua triangles), 
50% (pink triangles), and 100% (purple triangles, complete PBC) randomly selected cross-boundary links are applied, (b), Ejfect of communication 
wormhole on control: «ex versus v when a wormhole link is applied to any pair of nodes closer than d^w = 13-0, 12.5, 12.3, and 12.0 for p^w = 0.4%, 2%, 
3%, and 5%, respectively. 



To intuitively understand the non-monotonous phenomenon, we 
name the nodes within the boundary region (region between the 
boundary and the FEE) as boundary nodes (BNs), the nodes sitting 
on the FEE as frame nodes (FNs), and the nodes within the FEE as 
inside nodes (INs). Each node moves continuously in the space and 
acts as one of the three types of nodes from time to time. Due to the 
cutoff in the communication circles in the boundary region, the BNs 
have lower average degree than the FNs and INs, namely (/cbn) < 
C^in) ~ {^¥n)- At high velocity, most BNs are in fact originally INs 
from the inside carrying (w^^) packages, where (wu^) is the average 
flow of the INs. Apparently, we have {wij^)/{kij^) < {wi^)/{k^j^). 
Thus, for BNs, each of their links delivers more packages out than 
the INs. As recipients, FNs and INs have the same average degree 
((^in) = {^vn)) while FNs have more BN neighbors than INs have. 
We then have (wfn) > {'^m)> i-e-. FNs receive more packages. In 
addition, the FNs have more neighbors than the BNs, and thus we 
have (wpn) > (wbn)- This mechanism produces the FEE with higher 
density of extreme events due to a sudden increase in the flow when 
BNs and INs move into the region of FNs. However, the FEE so 
generated also increases the threshold 1^,3 and accordingly decreases 
the probability of occurrence of extreme events. Consequently, the 
coexistence of the two competing mechanisms generates the non- 
monotonous behavior, leading to an optimal velocity. Our theory, 
albeit non-rigorous, quantifies the relative strength of the two 
mechanisms and determines whether extreme events are facilitated 
or suppressed as velocity is increased. 

The effect of wormholes. In fact, consider an artificial system in 
which a fraction of links can be made cross-boundary. We find 
that, as the number of such "periodic-boundary-condition type" of 
links is increased, the ability for the mobfle network to suppress 
extreme events is degraded dramatically, as shown in Fig. 4a. This 
finding inspires significant applications in the emerging field of 
cybersecurity. Here we present one example by arguing that our 
control strategy can be used to detect wormhole attacks on mobile 
networks. In particular, a communication wormhole is a long-range 
link connecting two nodes beyond each's communication range, a 
recently developed concept in the field of network security''' 
Suppose any pair of nodes with distance over are connected by 
a wormhole. A smaller value of c/„h wUl then introduce more 
wormholes into the system. Let p„h be the ratio between the 
number of wormholes and the total number of edges E in the 
network. As shown in Fig. 4b, for moderate or relatively high 
velocity, a slight decrease in d^^ can trigger an explosive increase 



in «ex but there is only incremental change in p„h- Thus, abnormally 
large number of extreme events in a mobile network are strong 
indication that the network is under wormhole attack. Likewise, 
instilling a small number of wormholes into an adversary mobile 
network can effectively destroy its functions through the 
generation of a large number extreme events. 

Suppression of extreme events in intermittently mobUe networks. 

A network can be made mobile but only intermittently by 
introducing a "frozen" time interval AT during which the network 
structure is fixed. Specifically, each time after certain nodes in the 
network have moved, the network structure remains invariant for AT 
time steps before the next movement. For AT = 0, the network 
structure changes at every time step. Regardless of whether nodes 
are moving, the packet flow on the network is not interrupted. As 
shown in Figs. 5a-d, the number n^^ of extreme events exhibits a 
non- monotonous behavior with velocity v (for fixed values of AT) or 
with AT (for fixed value of v). In particular. Figs. 5a,b show, under 
hard-waU boundary conditions, «ex versus v for a number of fixed 
values of AT and versus AT for three values of the velocity, 
respectively. We observe that, increasing AT slightly above 0 leads 
to a sharp decrease in n^^, indicating that making the network 
intermittenfly mobile can be a cost-effective strategy to suppress 
extreme events. This holds regardless of the boundary condition. 
In fact, an intermittently mobile network system under periodic 
boundary conditions possesses a similar ability to suppress 
extreme events, as shown in the corresponding plots in Figs. 5c,d. 
All these indicate mobility as a general and effective strategy to 
suppress extreme events. 

While the concept of "mobility" can readily be generalized to static 
networks in which the nodes are not allowed to move physically, the 
existence of a physical boundary plays a critical role in controlling 
extreme events. For example, a network can be regarded as "mobile" 
insofar as its connection topology is time-varying, which can be 
realized, for example, by random rewiring of a subset of links from 
time to time in a computer network. In this case, an actual physical 
domain for the network cannot be properly defined. To demonstrate 
the importance of a physical boundary, we consider a classic ER 
network^^ with random rewiring. Let be the rewiring probabflity, 
the probability that an arbitrary link in the network is removed while 
a new link is established between a randomly selected pair of nodes 
that were not previously linked. Since the ER random network has no 
spatial structure, it is the change of nodal degree that leads to the 
variation in n^^ as is increased. As shown in Fig. 5e, Hex cannot be 
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Figure 5 | Number of extreme events in intermittently mobUe networks, (a) and (b), For a mobile network with hard-wall boundary conditions, ?3ex 
versus v and for a number of fixed values of the frozen time AT and versus AT for three values of v. (c) and (d), Corresponding results for periodic 
boundary conditions, (e) and (f), For a classic ER-network under intermittent random rewiring (a kind of generalized mobility), n^^ versus the rewiring 
probability Pj; (for a number of fixed values of AT) and versus AT (for Pjt = 0.5, 0.8, and 1.0), respectively. 



minimized by Pjj, indicating inability for random rewiring to sup- 
press extreme events. Simulations also show that the number of 
extreme events is not affected by the detailed rewiring mechanism, 
insofar as the nodal degree is unchanged. However, intermittent 
rewiring can reduce Wex to certain extent, as shown in Fig. 5f. For 
high value of in the ER random network, the number of extreme 
events, n^, can be solved analytically via the approach of self- 
consistent equation (see Methods Section). 



Degree bias and scaling. A previous result^" indicated that, in static 
networks, extreme events tend to occur more frequently on low- 
degree nodes. However, in any real-world network, certain degree 
of mobility is expected. Our computation indicates that, insofar as 
the network is mobile, extreme events tend to occur on highly 
connected nodes whose degrees are larger than the average degree 
of the network, as shown in Fig. 6a for a mobile network with hard- 
wail boundary conditions. Similar behaviors have been observed for 




80100 



Figure 6 | Degree bias and scaling. For a mobile network in a square domain with hard-wall boundary conditions, (a), number n^^ of extreme events 
occurring on small (red circles) and large (blue squares) degree nodes versus v, and (b), versus the total number of packets Wfor various velocity 
values. 
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mobile networks with periodic boundary conditions and for ER- 
random networks as well. We also find a scaling phenomenon 
between n^^ and the total number of packets W, as exemplified in 
Fig. 6b. 

Discussion 

To summarize, we find that extreme events in a complex network can 
be effectively controlled by making the network mobile, where the 
concept of mobility is general and can be implemented in various 
ways. A straightforward manifestation is that every node moves at a 
certain velocity but in random directions in a finite physical domain, 
in which case an analytic theory is developed to understand the 
existence of an optimal velocity to minimize the number of extreme 
events. Variants of the mobile schemes can also be considered, prov- 
ing the generality of our control strategy. Our theory has significant 
applications in the field of cybersecurity, where we articulate and 
demonstrate, for example, that wormhole attack on a mobile network 
can be detected through the monitoring of extreme events. Control of 
extreme events is an extremely important and challenging problem 
in science and engineering, and our work may stimulate further 
efforts. 

Methods 

Relation between ^(v) and T((v). Figure 7 show the relationship between a node's flow 
w{t — 1) at time t— I and its degree at time f for different velocity values, the basic 
relation employed in our theoretical analysis in the main text, and the function ?/(v) is 
shown in Fig. 7d. We observe a robust linear relation for mobUe networks with 
periodic boundary conditions and for ER-random networks^^. For mobile networks 



with hard-wall boundary conditions, the relation is still linear, especially for low or 
high velocity values. For moderate velocity values, there is deviation from the linear 
behavior, due to the boundary effect as shown in Fig. 3. {This boundary effect does not 
exist in mobile networks with periodic boundary conditions and in ER-random 
networks.) In general, we observe an approximately linear relation between w{t — 1) 
and k{t), which can be taken as an empirical law for our analysis of the extreme-event 
dynamics in mobile networks. 

We present a detailed derivation of the relation between (i{v) and //{v). The starting 
point is the following packet conservation law: 



W = N ^ {w{t-\))P(k) 



W 



1(0 = 0 

N-l 



= E m-k{t)+n{v)]p{k{t)) 



(3) 



k{l)=0 



t(<)=o 



W 



- = c{vym)+,^(v) 



The average degree of the whole network is {k{t)) — 2E/N. We thus have 

W/N-ij{v) _ W-N-ri{v) 
2E/N ~ 2E 



cM = - 



(4) 



Detailed analytical explanation to the non-monotonous behavior. As illustrated in 
Fig. 2a, the probability that z packets are delivered to the focal node located at b from a 
neighbor of degree k^it) located at and carrying w^it — 1) packets is given by 
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Figure 7 | Linear relation between {w(f— l))andic(f). (a-c). Linear relation between (w(t— 1)) and k{t) for mobile networks with hard-wall and periodic 
boundary conditions, and for ER-random networks, respectively, (d). Dependence of );(v) on v for mobile networks with hard-wall and periodic 
boundary conditions (PBCs), and )/(Pk) versus Pr for ER-random networks (inlets). All results are averaged over 20 statistical realizations. 
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P{z,K{t),w„{t-l)) = 



w„{t-l) 



1 



(5) 



■p,{K{t)yp{w,{t-i)), 



where Pkiknit)) is the degree distribution of the neighboring nodes: 

1 N-2-k„{t) 



PkiUt))-- 



N-2 
knit) 



12 



and S{bn) is the circular communication area of the neighbor located at b^. For ^ LI 
2 — a, the area is a complete circle with S(b^) — na^. However, forL/2 — a<b,^^ 1/2, 
a cutoff by the boundary leads to 5(feiJ — 71^^(1 — a/n) + asina ■ [Lll — where tx 
= arccos{(L/2 — b^^ld), as shown in Fig. 2a. The quantity Piw^it — 1)) is the 
probability density function (PDF) of wj^t — 1), which depends on the velocity v and 
further analysis is needed to determine its mathematical form. 

Using Eq. (1), we can write the probability for a packet to visit a node as {w)!W — 
(1/ WQ [i;(0) ■ /c + );(0)] for v = 0. For v # 0, the probability for a node at fj„ to be visited 
isp„ = {wjyt - l))IW = [c(v) ■ kj^t) + ri{v)]IW. We thus have 

Substituting this expression into Eq. (5), we obtain the PDF of z as 

V — \W fv — \N — 2 ~| 

^(^) = Z^«.„(,-i)=z [Z^t.(,) = i Ut),Mt-'i))\ ■ Consider any one of the 
boundaries of the square domain, as shown in Fig. 2a. For the focal node located at b, 
Z = - i^i pa^ckets will be delivered to it from its m neighbors located in the bar- 
shaped region of width Ab^^ ^ 0 at b^. The PDF of Z associated with is the 
convolution of the m PDFs, taking into account all possible values of m: 

PbA^) = Ylm-Q W^i) * -^^^2) * ... * P{Zni)]-P(m), where P(m) is the probability 
for m nodes to exist in the bar region: 



PK(t-i))=El! 



P(m) 



JV-1 



"17 



(N-l-m) 



(6) 



Abn- The PDF of w{t) can then be expressed as the 



with area given by Sbar — 2Aba ^J~^^--{ba--bf' for small A^n- The flow w{t) of the 
node at time t is the sum of the incoming packets from nodes in all the bar-shaped 
regions within its communication circle: w(f) — Zj, where the contribution 

from each region can be treated independently and 
|minQ, b + a^ -{h- 

convolution of all Pj,^ (■2j)'s for e [b — a, min{L/2, b + a)], as 
Tfe(w(f)) ^Pb„, (Zi) i Pb„, * • • • * P,„, (Zm). 

For a static node located at bo, the PDF of its flow is Tj,^ {w{t)). For any moving 
node (v > 0) initially at bo, the PDF of its flow is contributed by the *Pf,(w} fields at all 
the locations it has visited, i.e., the circular region centered at bo with estimated radius 
R{v). The PDF of this node's flow thus contains the contributions from all Tt(w(t))'s 
for b ^ [bo — Riv), min{fco + R{v),L/2)]. The random walk radius i^{v) increases with 
V, andi^{v) Ofor 0 and reaches L/2 {the upper bound of /?) for v^'~-^. Since 
and rjiv) both increase monotonically with v, we have, approximately, 
R(v) c^rj{v)/ {2W/NL). Consequently, the flow PDF of a node initially at location b 
moving with velocity v can be obtained as a function of b: 



pmm[b+R(v),i/2] 
Jmax[0,b-R{v)] 



4'i,(w(0)-^dfo 



|-mlnl6+Ji(v),L/2| 
Jmax[0,fj-R(v; 



(7) 



rdb 



where 8b dblV' is the probability that a node is in the region about the location b via 
random walk, and the denominator on the right side is a normalization factor. The 
quantity ^^{w{t)) determines the threshold — {w{t)),-^ + ^^w^ over which 

V — \W 

extreme events take place, where {w{t))^^^ — MO'^'b(''^(0) ^^"^ 

5^ — [MO ~ (MO)'-^]'^- Finally, we obtain the number of extreme 

events occurred within T time steps as 



.L/2 

Jo 



«'(0=<!„ + l 



(8) 



which is a function of ;^(v}. 

In the derivation of Eq. (2), the degrees of the neighboring nodes are treated 
independently in Eq. (5) and the expression of P(w„(l — 1)), which is an approxi- 
mation because nodes that are spatially close tend to be degree-correlated. The spatial 
degree correlation represents the main discrepancy between our moving random- 
network model and a mobile network in the real physical space (see the following for a 



detailed discussion about the effect of spatial correlation). Notwithstanding, the 
degree correlation has no dependence on the velocity and, hence, there is little effect 
on extreme-event generation in the system for different values of the velocity. 

The average flow {w{t))[, at location b can be obtained via T(,(w(i)) as well as the 
mean -field theory: 



Jk—a 



L2 



dbr, 



(Kit)) 

na'-[\->^+a{\-a-b)sml^ {w,(t-\)) 
{K{t)) 



(9) 



L2 



where ^ = arccos[(L/2 - a- h)la\, {wj,t - 1)> = ^{v) - {k^{t)) + ?/(v), and 
{kn{t)) — k„{t)'Pk{kn{t)). The first term on the right side is the contri- 

bution from the neighboring nodes with their communication circles partially cutting 
off by the boundary, while the second term represents the contribution from neigh- 
bors with complete communication circles. 

Self-consistent equation to estimate the number of extreme events in ER random 
networks. For high value of P^ in the ER random network, the number of extreme 
events, Mex> can be solved analytically via the approach of self-consistent equation. 
Regarding P^ as corresponding to some kind of "velocity" v, we have that the number 
of extreme events n^^ for v — 0 is given by 



E 

^=1+1 



w 



IE) \ IE 



Pk{k), 



(10) 



where q is the threshold by which extreme events are defined, and the degree 
distribution P{k) of the network can be written in a binomial form: 



N-1 



(11) 



To map the ER random network into a two-dimensional mobile network, we set the 
link probability to be p — na^lV. Eq. (lO)yields «ex — 125, which matches the 
simulation result for v — 0 well. 

For high velocity, the PDF of flow P{w) becomes independent of k and identical for 
all nodes, which can be solved analytically by the self-consistent equation method. In 
particular, the number z of packets delivered to a node of degree k from one of its 
neighbors of degree fc^ has the following PDF 



f'W=E 



E^(" 



Pkikn), 



(12) 



where P^(/:n) is the degree distribution of the neighboring node given by Eq. (11), and 
P{w) can be obtained via the weighted sum of the convolution of P{z): 



K^) - I-^t^i) * -^(^2) * ■ ■ ■ * P{zk)]-Pk{k). 



(13) 



V — \W 

By substituting any form of P{w) that satisfies _q w-P{w) — W/N into Eq. (13), 
we obtain an updated form of P{w) via Eq. (12). The solution of P(w) can be obtained 
by iterating the self-consistent equation until the root- mean -square deviation 
between P{w) and its updated version falls below a small threshold e (e.g., €— 10~*). 
Figure 8a shows a good agreement between P(w) from direct simulation and from the 
self-consistent approach. The number of extreme events can then be estimated, 



yielding "ex^NT^^^^^ 
direct simulation. 



P(w) ^ 238, which agrees quite well with the value from 



Effect of degree-degree correlation. The probability pA that a given node i and one of 
its neighbors j share a neighbor in the ER- random network of link probability p isp^. 
For a mobile network in a two-dimensional domain, this probability is given by 



PA- 



r Sij{r) 27ir 
Jo na^ L2 



dr, 



(14) 



where r denotes the distance between node / and its neighboring node j, and Sjj (r) is 
the overlapping area of /' and/s communication circles, which is given by 



(15) 



Sij{r)^a arccos(^— j-rya2__ 



Numerical calculation reveals thatpA^-p , indicating that a node's degree is much 
more likely to have a value similar to that of its neighbor, giving rise to a strong 
correlation between the degrees of neighboring nodes, or the phenomenon of degree- 
degree correlation. Namely, node /"s degree k is correlated with its neighbor's degree 
ky^, the so-called k-k^ correlation. In addition, the degrees of the neighboring nodes of i 
are also correlated, leading to the so-called k^-k^ correlation. Since the effect of cutoff 
of the communication circle by a boundary makes the analysis on spatial degree 
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Figure 8 | PDF of flow P{w) in simulation and in different theory assumptions, (a), For the classic ER random networks with rewiring probability Pj^ — 
1, PDF of flow P{w) obtained from direct simulation (red squares) and self-consistent equation (solid blue line), (b), PDF of flow P{w) from direct 
simulation (red squares), from theory without taking into account spatial degree correlation (black circles), and from theory with degree-degree 
correlation (solid blue line). 



correlation extraordinarily complicated, liere we neglect the effect of the boundary or, 
equivalently, we focus on mobile networks with periodic boundary conditions. 

For convenience, we classify the neighbors of node j into two classes with respect to 
the focal node i: riji neighbors in the overlapping area Sy(r) (i.e., the regions shared by 



these nodes and /) and the remaining 



area. The PDF of n,i is 



■ Hji neighbors in the non- overlapping 



(16) 



wherepi — 5,-, {r)l{na^). The number tt^2of neighbors in the non- overlapping area has 
the PDF as 



N-ki{l)-2-n 



where p2 ~ {t^i^^ ~ ^iji^))!^^- The degree distribution of node j is given by 

Inr . 



Jo -TTT- 



-dr 



(17) 



(18) 



";l+n>2=fcn 



where /c, — n^y -I- np_. In the high-velocity regime, the self-consistent equation method 
can be adopted to solving P(w) by substituting Pkikn) from Eq. (12) into Eq. (18), 
while the quantity Pk{k) in Eq. (13) can still be determined by Eq. (11). Figure 8b 
shows the flow distribution P{w) obtained from direct simulation and from self- 
consistent theory with or without the k-k^ correlation. We observe a better agreement 
between simulation and theory when the k-k^ correlation is taken into account. 

It is mathematically difficult to consider both the k-k^^ and /Cn-^n correlations. 
Nonetheless, since the spatial degree correlation depends solely on the spatial dis- 
tribution of the nodes, which is completely random regardless of the velocity, the 
combined effect of the two kinds of correlations on has little dependence on the 
actual value of the velocity. Consequently, the non-monotonous behavior of extreme 
events caused by the boundary effect should persist even in realistic situations where 
both types of degree-degree correlations are present, and we expect the theory of 
mobility-based control developed in the main text to be generally valid. 
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